The fidelity of dynamic signaling by noisy biomolecular networks

This is the final version of the article. Available from Public Library of Science via the DOI in this record.Cells live in changing, dynamic environments. To understand cellular decision-making, we must therefore understand how fluctuating inputs are processed by noisy biomolecular networks. Here we present a general methodology for analyzing the fidelity with which different statistics of a fluctuating input are represented, or encoded, in the output of a signaling system over time. We identify two orthogonal sources of error that corrupt perfect representation of the signal: dynamical error, which occurs when the network responds on average to other features of the input trajectory as well as to the signal of interest, and mechanistic error, which occurs because biochemical reactions comprising the signaling mechanism are stochastic. Trade-offs between these two errors can determine the system's fidelity. By developing mathematical approaches to derive dynamics conditional on input trajectories we can show, for example, that increased biochemical noise (mechanistic error) can improve fidelity and that both negative and positive feedback degrade fidelity, for standard models of genetic autoregulation. For a group of cells, the fidelity of the collective output exceeds that of an individual cell and negative feedback then typically becomes beneficial. We can also predict the dynamic signal for which a given system has highest fidelity and, conversely, how to modify the network design to maximize fidelity for a given dynamic signal. Our approach is general, has applications to both systems and synthetic biology, and will help underpin studies of cellular behavior in natural, dynamic environments.We acknowledge support from a Medical Research Council and Engineering and Physical Sciences Council funded Fellowship in Biomedical Informatics (CGB) and a Scottish Universities Life Sciences Alliance chair in Systems Biology (PSS). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript

Histone Modifications at Human Enhancers Reflect Global Cell-Type-Specific Gene Expression

The human body is composed of diverse cell types with distinct functions. Although it is known that lineage specification depends on cell-specific gene expression, which in turn is driven by promoters, enhancers, insulators and other cis-regulatory DNA sequences for each gene1, 2, 3, the relative roles of these regulatory elements in this process are not clear. We have previously developed a chromatin-immunoprecipitation-based microarray method (ChIP-chip) to locate promoters, enhancers and insulators in the human genome4, 5, 6. Here we use the same approach to identify these elements in multiple cell types and investigate their roles in cell-type-specific gene expression. We observed that the chromatin state at promoters and CTCF-binding at insulators is largely invariant across diverse cell types. In contrast, enhancers are marked with highly cell-type-specific histone modification patterns, strongly correlate to cell-type-specific gene expression programs on a global scale, and are functionally active in a cell-type-specific manner. Our results define over 55,000 potential transcriptional enhancers in the human genome, significantly expanding the current catalogue of human enhancers and highlighting the role of these elements in cell-type-specific gene expression

Transcription factor site dependencies in human, mouse and rat genomes

<p>Abstract</p> <p>Background</p> <p>It is known that transcription factors frequently act together to regulate gene expression in eukaryotes. In this paper we describe a computational analysis of transcription factor site dependencies in human, mouse and rat genomes.</p> <p>Results</p> <p>Our approach for quantifying tendencies of transcription factor binding sites to co-occur is based on a binding site scoring function which incorporates dependencies between positions, the use of information about the structural class of each transcription factor (major/minor groove binder), and also considered the possible implications of varying GC content of the sequences. Significant tendencies (dependencies) have been detected by non-parametric statistical methodology (permutation tests). Evaluation of obtained results has been performed in several ways: reports from literature (many of the significant dependencies between transcription factors have previously been confirmed experimentally); dependencies between transcription factors are not biased due to similarities in their DNA-binding sites; the number of dependent transcription factors that belong to the same functional and structural class is significantly higher than would be expected by chance; supporting evidence from GO clustering of targeting genes. Based on dependencies between two transcription factor binding sites (second-order dependencies), it is possible to construct higher-order dependencies (networks). Moreover results about transcription factor binding sites dependencies can be used for prediction of groups of dependent transcription factors on a given promoter sequence. Our results, as well as a scanning tool for predicting groups of dependent transcription factors binding sites are available on the Internet.</p> <p>Conclusion</p> <p>We show that the computational analysis of transcription factor site dependencies is a valuable complement to experimental approaches for discovering transcription regulatory interactions and networks. Scanning promoter sequences with dependent groups of transcription factor binding sites improve the quality of transcription factor predictions.</p

Aurora-A Interacts with AP-2α and Down Regulates Its Transcription Activity

Aurora-A is a serine/threonine protein kinase and plays an important role in the control of mitotic progression. Dysregulated expression of Aurora-A impairs centrosome separation and maturation, which lead to disrupted cell cycle progression and tumorigenesis. However, the molecular mechanism by which Aurora-A causes cell malignant transformation remains to be further defined. In this report, using transcription factors array and mRNA expression profiling array, we found that overexpression of Aurora-A suppressed transcription activity of AP-2α, a tumor suppressor that is often downregulated in variety of tumors, and inhibited expression of AP-2α-regulated downstream genes. These array-based observations were further confirmed by microwell colorimetric TF assay and luciferase reporter assay. Downregulated transcription activity of AP-2α by Aurora-A was found to be associated with reduced AP-2α protein stability, which appeared to be mediated by Aurora-A enhanced ubiquitin-dependent proteasomal degradation of AP-2α protein. Interestingly, Aurora-A-mediated AP-2α degradation was likely dependent Aurora-A kinase activity since inhibition of Aurora-A kinase activity was able to rescue Aurora-A-induced degradation of AP-2α. Moreover, we defined a physical interaction between Aurora-A and AP-2α, and such interaction might bridge the suppressive effect of Aurora-A on AP-2α protein stability. These findings provide new insights into molecular mechanism by which Aurora-A acts as an oncogenic molecule in tumor occurrence and malignant development

Trade-Offs and Constraints in Allosteric Sensing

Sensing extracellular changes initiates signal transduction and is the first stage of cellular decision-making. Yet relatively little is known about why one form of sensing biochemistry has been selected over another. To gain insight into this question, we studied the sensing characteristics of one of the biochemically simplest of sensors: the allosteric transcription factor. Such proteins, common in microbes, directly transduce the detection of a sensed molecule to changes in gene regulation. Using the Monod-Wyman-Changeux model, we determined six sensing characteristics – the dynamic range, the Hill number, the intrinsic noise, the information transfer capacity, the static gain, and the mean response time – as a function of the biochemical parameters of individual sensors and of the number of sensors. We found that specifying one characteristic strongly constrains others. For example, a high dynamic range implies a high Hill number and a high capacity, and vice versa. Perhaps surprisingly, these constraints are so strong that most of the space of characteristics is inaccessible given biophysically plausible ranges of parameter values. Within our approximations, we can calculate the probability distribution of the numbers of input molecules that maximizes information transfer and show that a population of one hundred allosteric transcription factors can in principle distinguish between more than four bands of input concentrations. Our results imply that allosteric sensors are unlikely to have been selected for high performance in one sensing characteristic but for a compromise in the performance of many

Adaptation-Dependent Synchronous Activity Contributes to Receptive Field Size Change of Bullfrog Retinal Ganglion Cell

Nearby retinal ganglion cells of similar functional subtype have a tendency to discharge spikes in synchrony. The synchronized activity is involved in encoding some aspects of visual input. On the other hand, neurons always continuously adjust their activities in adaptation to some features of visual stimulation, including mean ambient light, contrast level, etc. Previous studies on adaptation were primarily focused on single neuronal activity, however, it is also intriguing to investigate the adaptation process in population neuronal activities. In the present study, by using multi-electrode recording system, we simultaneously recorded spike discharges from a group of dimming detectors (OFF-sustained type ganglion cells) in bullfrog retina. The changes in receptive field properties and synchronization strength during contrast adaptation were analyzed. It was found that, when perfused using normal Ringer's solution, single neuronal receptive field size was reduced during contrast adaptation, which was accompanied by weakening in synchronization strength between adjacent neurons' activities. When dopamine (1 µM) was applied, the adaptation-related receptive field area shrinkage and synchronization weakening were both eliminated. The activation of D1 receptor was involved in the adaptation-related modulation of synchronization and receptive field. Our results thus suggest that the size of single neuron's receptive field is positively related to the strength of its synchronized activity with its neighboring neurons, and the dopaminergic pathway is responsible for the modulation of receptive field property and synchronous activity of the ganglion cells during the adaptation process

Large Scale Gene Expression Profiles of Regenerating Inner Ear Sensory Epithelia

Loss of inner ear sensory hair cells (HC) is a leading cause of human hearing loss and balance disorders. Unlike mammals, many lower vertebrates can regenerate these cells. We used cross-species microarrays to examine this process in the avian inner ear. Specifically, changes in expression of over 1700 transcription factor (TF) genes were investigated in hair cells of auditory and vestibular organs following treatment with two different damaging agents and regeneration in vitro. Multiple components of seven distinct known signaling pathways were clearly identifiable: TGFβ, PAX, NOTCH, WNT, NFKappaB, INSULIN/IGF1 and AP1. Numerous components of apoptotic and cell cycle control pathways were differentially expressed, including p27KIP and TFs that regulate its expression. A comparison of expression trends across tissues and treatments revealed identical patterns of expression that occurred at identical times during regenerative proliferation. Network analysis of the patterns of gene expression in this large dataset also revealed the additional presence of many components (and possible network interactions) of estrogen receptor signaling, circadian rhythm genes and parts of the polycomb complex (among others). Equal numbers of differentially expressed genes were identified that have not yet been placed into any known pathway. Specific time points and tissues also exhibited interesting differences: For example, 45 zinc finger genes were specifically up-regulated at later stages of cochlear regeneration. These results are the first of their kind and should provide the starting point for more detailed investigations of the role of these many pathways in HC recovery, and for a description of their possible interactions

A functional SNP in the regulatory region of the decay-accelerating factor gene associates with extraocular muscle pareses in myasthenia gravis

Complement activation in myasthenia gravis (MG) may damage muscle endplate and complement regulatory proteins such as decay-accelerating factor (DAF) or CD55 may be protective. We hypothesize that the increased prevalence of severe extraocular muscle (EOM) dysfunction among African MG subjects reported earlier may result from altered DAF expression. To test this hypothesis, we screened the DAF gene sequences relevant to the classical complement pathway and found an association between myasthenics with EOM paresis and the DAF regulatory region c.-198C>G SNP (odds ratio=8.6; P=0.0003). This single nucleotide polymorphism (SNP) results in a twofold activation of a DAF 5′-flanking region luciferase reporter transfected into three different cell lines. Direct matching of the surrounding SNP sequence within the DAF regulatory region with the known transcription factor-binding sites suggests a loss of an Sp1-binding site. This was supported by the observation that the c.-198C>G SNP did not show the normal lipopolysaccharide-induced DAF transcriptional upregulation in lymphoblasts from four patients. Our findings suggest that at critical periods during autoimmune MG, this SNP may result in inadequate DAF upregulation with consequent complement-mediated EOM damage. Susceptible individuals may benefit from anti-complement therapy in addition to immunosuppression

Evolution of protein domain architectures

This chapter reviews current research on how protein domain architectures evolve. We begin by summarizing work on the phylogenetic distribution of proteins, as this will directly impact which domain architectures can be formed in different species. Studies relating domain family size to occurrence have shown that they generally follow power law distributions, both within genomes and larger evolutionary groups. These findings were subsequently extended to multi-domain architectures. Genome evolution models that have been suggested to explain the shape of these distributions are reviewed, as well as evidence for selective pressure to expand certain domain families more than others. Each domain has an intrinsic combinatorial propensity, and the effects of this have been studied using measures of domain versatility or promiscuity. Next, we study the principles of protein domain architecture evolution and how these have been inferred from distributions of extant domain arrangements. Following this, we review inferences of ancestral domain architecture and the conclusions concerning domain architecture evolution mechanisms that can be drawn from these. Finally, we examine whether all known cases of a given domain architecture can be assumed to have a single common origin (monophyly) or have evolved convergently (polyphyly). We end by a discussion of some available tools for computational analysis or exploitation of protein domain architectures and their evolution